The nature of Feshbach molecules in Bose-Einstein condensates 
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We discuss the long range nature of the molecules produced in recent experiments on molecular 
Bose-Einstein condensation. The properties of these molecules depend on the full two-body Hamil- 
tonian and not just on the states of the system in the absence of interchannel couplings. The very 
long range nature of the state is crucial to the efficiency of production in the experiments. Our 
many-body treatment of the gas accounts for the full binary physics and describes properly how 
these molecular condensates can be directly probed. 
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Bose-Einstein condensation of molecules is an excit- 
ing challenge and opportunity in the physics of ultracold 
gases. As direct laser cooling of molecules is largely pre- 
vented by their densely lying rovibrational energy lev- 
els, several approaches now focus on the association of 
atoms in Bose-Einstein condensates. Present techniques 
are based on photoassociation Q] or magnetic field tun- 
able interactions 0, Q . The detection of molecular con- 
densates, however, remains a demanding problem. Re- 
cent experiments at JILA Q found evidence for molecu- 
lar condensation by probing coherence properties of the 
assembly of atoms plus molecules. The 85 Rb conden- 
sate was exposed to a sequence of two fast magnetic field 
pulses in the vicinity of a Feshbach resonance. The pulses 
were separated by a variable period of time with a sta- 
tionary magnetic field, termed the evolution period. The 
observed final densities of gas atoms, i.e. the remnant 
condensate and a "burst" of comparatively hot atoms, in- 
dicated a coherent coupling between atoms and diatomic 
molecules in a highly excited vibrational state during the 
evolution period. Several subsequent theoretical stud- 
ies 0, El have concluded that the gas contained a 
molecular condensate at the end of the pulse sequence. 
The magnitude of the molecular fraction, and the pre- 
cise mechanism of producing the molecules are a mat- 
ter of continuing controversy d, . To settle this issue 
the next generation of experiments could determine the 
coherent superposition of atomic and molecular compo- 
nents by directly detecting the molecules. The efficiency 
of different detection schemes based on laser excitation 
|ic| depends sensitively on the number of molecules and 
the wave functions of the diatomic bound states produced 
by a magnetic field pulse. 

In this letter we shall give a full description of these 
wave functions for typical experimental magnetic field 
strengths and determine the molecular component in the 
evolution period of the pulse sequence. In particular, 
we shall show explicitly that the wave functions of the 
molecules have a spatial extent of the order of the scat- 



tering length. At the magnetic field strengths closest 
to resonance this length scale even becomes comparable 
with the mean distance between the atoms in the di- 
lute condensate. Under these conditions a separation of 
the gas into atoms and diatomic molecules is physically 
meaningless. The strong binary correlations provided by 
the long range intermediate molecular states, however, 
are still crucial for the efficiency of the association of 
atoms to molecules. In the evolution period the typi- 
cal static experimental field strengths are sufficiently far 
from the resonant field that the highly excited diatomic 
molecules can exist as a metastable entity of the gas. The 
microscopic many-body appr oach we apply to determine 
the molecular fraction accounts for the long range 

nature of the intermediate molecular states. This ap- 
proach treats bound and free molecular states formed 
during the pulse sequence in a unified manner and pro- 
vides a straightforward understanding of the association 
of condensate atoms to molecules. 

The experimental technique of magnetic field tunable 
interactions takes advantage of the Zeeman effect in the 
electronic energy levels of the atoms. In the JILA ex- 
periments 0| the 85 Rb condensate atoms were prepared 
in the (F = 2, nip — —2) hyperfine state. Throughout 
this letter the open s-wave binary scattering channel of 
two asymptotically free atoms in the (F = 2, mp = —2) 
state will be denoted as the {—2, —2} open channel with 
an associated reference potential Vh g (r). When the mp 
degeneracy of the hyperfine levels is removed by an ex- 
ternal magnetic field B the potentials associated with 
the different asymptotic scattering channels are shifted 
with respect to each other. Although the {—2, —2} open 
channel is only very weakly coupled to other open chan- 
nels, it can be strongly coupled to closed channels. A 
zero-energy scattering resonance occurs when the field- 
dependent energy E rcs (B) of a closed channel vibrational 
state (a Feshbach resonance level) <^ res is tuned close to 
the dissociation threshold energy of Vb g (r). We note that 
the closed channel state <^ res is not a stationary state of 
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FIG. 1: Scheme of a typical magnetic field pulse shape in the 
low density (no = 3.9 x 10 12 cm -3 ) experiments in Ref. |4(- 
The minimum magnetic field strength of the first and second 
pulse is -B m in = 15.55 mT. In the evolution period the field 
strength is chosen as -Bevoivc = 16.0 mT. In the course of the 
experiments the evolution time t cvo \ vc as well as -Bevoive were 
varied. The dashed line indicates the position of the resonance 
at B = 15.49 mT. 



the full two-body Hamiltonian, which has a coupling be- 
tween the channels. If E ICS (B) approaches the threshold 
from below, however, the overall potential matrix sup- 
ports a shallow multi-channel bound state 0b with energy 
i?b- This proper stationary molecular state ceases to ex- 
ist at the position of the resonance {Bq — 15.49 mT Q), 
for which the s wave scattering length a of two asymp- 
totically free atoms in the {—2, —2} open channel has 
a singularity. When E\, is sufficiently close to thresh- 
old, the component of 0b in the 1—2, —2} open channel 
then exhibits the universal form [12j exp(— r/a)/r at large 
relative distances of the two atomic constituents of the 
molecular state. 

We shall now describe more fully the long range nature 
of the molecular states for typical experimental magnetic 
field strengths in the fast sequence of pulses shown in 
Fig. ^ To this end we have reduced the complete multi- 
channel Hamiltonian of the relative motion of two atoms 
to a two-channel model: 



H 



2B 



W(r) 



V hs (r) W(r) 

-^V 2 + V cl (B,r) 



(1) 



Here m is the atomic mass of 85 Rb. For the potential in 
the {—2, —2} open channel we use a Lennard-Jones form 
Vbg(r) = 4e [(cr/r) 12 - (a/r) 6 ] with a = 37.3292 aBohr 
(a Bo hr = 0.052918 nm) and 4ecr 6 = C 6 = 4660 a.u. [l3| 
(1 a.u. = 0.095734 yj nm 6 ). Vb g {r) then reproduces the 
background scattering length of Ob g = —450 aBohr 0- I n 
accordance with Ref. 0] we model the closed channel 
potential as V c \(B,r) = Vb g (r) + E C \(B), where E C \(B) 
follows the dependence of the energy difference of the 
corresponding Zeeman hyperfine levels on the magnetic 



field. We use h^dE^/dB = -34.6 MHz/mT. In this 
simplified model V c \{B,r) supports only two vibrational 
states. We assume the excited one to be the resonance 
state res which thus satisfies [— S, 2 V 2 /m + Vy</> rcs ~ 
Rgs^res- We have chosen the off diagonal potential as 
[3 W(r) = /3exp(-r/a) with f3/k B = 38.5 mK and 
a = 5 dBohr, which gives a width of the resonance of 
AB = 1.1 mT |4j. With this choice of parameters the 
model also produces the correct shift between the reso- 
nance position Bo and the magnetic field strength B ves , 
where the resonance state crosses the dissociation thresh- 
old of Vb g {r) (E les (B les ) = 0). For the 85 Rb resonance 
at B = 15.49 mT this shift is negative and of the or- 
der of Bq — B ICS = —0.9 mT, so that the multi-channel 
bound state 0b persists when the resonance state has 
crossed the threshold. The corresponding wave functions 
and negative binding energies are obtained from Eq. 
by -H2B</>b = -Eb^b and the scattering length is approxi- 
mated by 



a(B) = a hs 1 - 



AB 

B — Br 



(2) 



Figure[2]shows the bound state 0b at the magnetic field 
strengths of B C voivc = 16.0 mT and B min = 15.55 mT. 
While -Bevoive corresponds to the magnetic field during 
the evolution period of the pulse sequence in Fig. Q 
-S m in is the magnetic field strength closest to resonance. 
Both molecular states in Fig. |2 are dominated by their 
{—2, —2} open channel component which contributes 
95.3% (99.9%) of the total probability density at Bevoive 
(-B m i n ). For both magnetic field strengths the small 
closed channel component of 0b assumes the form of the 
resonance state 1OS (r) . Close to the resonance the bond 
length of 0b (mean internuclear distance) is on the or- 
der of (r) = a/2. One obtains (7-) = 326 aBohr (a = 
521 aBohr from Eq. J2J) at B cvo \ vc and (r) = 4255 asohr 
(a = 7800 aBohr from Eq. (0) at B m i n . At the lowest 
magnetic field strengths in the pulse sequence in Fig. ^ 
the spatial extent of the molecules thus becomes com- 
parable with the mean distance between the condensate 
atoms in the JILA experiments 0] , which is of the order 

— 1/3 

of n Q '" = 0.64 /im = 12000 aBohr- Under these con- 
ditions, the dilute gas parameter \/noa 3 is comparable 
to 1, and one can no longer identify a particular pair of 
atoms in the g molecule in the state 0b because 

its molecular wave function would overlap with other gas 
atoms. One and the same atom could thus contribute to 
several diatomic molecules. After the pulse sequence and 
during the evolution period, however, the gas is weakly 
interacting (noa 3 <^ 1) and diatomic molecular bound 
states are sufficiently confined in space for the number of 
molecules to be a meaningful quantity. 

We shall illustrate the resulting remarkable evolution 
of the molecular fraction in the JILA experiments with 
the quantum mechanical many-body description devel- 
oped in Ref. [jj for a homogeneous gas with the exper- 
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FIG. 2: Coupled-channel bound states at the magnetic field 
strengths of -B CV oive = 16.0 mT and B m i n = 15.55 mT. 
The radial coordinate is given on a logarithmic scale. The 
dotted (dashed) curves indicate the closed ({—2, —2} open) 
channel components. The solid curves are the correspond- 
ing bound state wave functions of the separable potential in 
Ref. 0- The separable potential wave functions agree with 
the {—2, —2} open channel components of the coupled chan- 
nels bound states at distances large compared to the van der 
Waals length of i(mC 6 /ft 2 ) 1/4 = 82 a Bohr . 

imental mean density of no = 3-9 x 10 12 cm~ 3 . The 
underlying cumulant approach 0, IT3 is based on the 
microscopic many-body Hamiltonian |jj and systemati- 
cally decouples the exact hierarchy of dynamic equations 
for quantum correlation functions at any desired degree 
of accuracy. According to Ref. the dynamics of the 
condensate mean field the molecular fraction of the 
gas, and the burst of correlated pairs of comparatively 
hot atoms 0, are all determined by a single nonlinear 
Schrodinger equation. In the homogeneous gas limit this 
has the form: 

f) r°° 8 

ih^*(t) = "**(*) J t dr ^( T )^-h(t,r). (3) 

Here to is the initial time of the pulse (see Fig. and 
hit, t) is the coupling function that accounts for the in- 
teractions of the gas atoms. In Ref. we have described 
the binary collision dynamics in terms of a single separa- 



ble potential V(t) that mimics the low energy scattering 
properties of the two component coupled-channel Hamil- 
tonian Q and, at each magnetic field strength B(t), the 
long range of the {—2,-2} component of the coupled- 
channel bound state (see Fig. [21 . The coupling function 
includes the entire two-body time development operator 
U2B (tj t) and is given by 

h(t,r) = (2nh) 3 0(t - T)(0\V(t)U 2B (t,T)\0). (4) 

Here 6(t—r) is the step function that yields 1 at t > r and 
vanishes elsewhere, and |0) is the zero momentum plane 
wave. In this treatment U2B(t,r) accounts, in a unified 
manner, for both bound and free molecular states formed 
during thepulse sequence. 

In Ref. |3 we have determined the number of diatomic 
molecules in a weakly interacting gas in terms of the 
quantum mechanical observable |l6( that counts, for each 
of the N(N — l)/2 pairs of atoms in a gas with N atoms, 
the overlap between the molecular bound state wave 
function and the many-body state. We have shown ex- 
plicitly in Ref. that, despite the N 2 scaling behavior of 
the number of pairs of atoms, our microscopic many-body 
description predicts the number of diatomic molecules to 
be smaller than N/2, as soon as noa 3 -C 1. According 
to Ref. 7], for a homogeneous gas, the amplitude of the 
density of the molecular fraction is given by 

*b(*) = 4 / d 3 r 0£(r) [$(r, t) + * 2 (t)] . (5) 

Here t is the time at which the molecules are detected 
(and thereby destroyed), and 

$(r,t) = -(27rft)i / dr ^ 2 ( T )^-(r\U 2B (t, r)|0) (6) 

J to 

is the pair function 0E1- The densities of the molecular 
fraction and the atomic condensate are then obtained by 
n b (f) = |* b (t)| 2 and n c (t) = |*(i)| 2 , respectively. 

Figure |3 shows the evolution of the atomic conden- 
sate density and the molecular fraction in a homogeneous 
gas with a density of no = 3.9 x 10 12 cm -3 for differ- 
ent evolution times of the magnetic field pulse sequence 
in Fig. n As explained above, the molecular fraction 
is physically meaningful only at the beginning (to = 0) 
and end (ta n = 84, 85, 86, 87 fxs in Fig. |3J of the pulse 
sequence as well as in the evolution period (beginning 
at 37 us in Fig. [SJ. At to = the gas is prepared as 
a pure atomic condensate so that nb(*o) = 0. The fi- 
nal densities n c (tfi n ) and n\y(tn n ) exhibit the experimen- 
tally observed Q pronounced oscillations depending on 
the evolution time t e voive (see Fig. , with a frequency 
\Eh{B cvo \ ve )\/h and a relative phase shift depending on 
the relative magnitudes of all components of the gas. The 
amplitude of the oscillation in nb(tfi n ) is about twice the 
near constant molecular fraction in the evolution period 
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FIG. 3: Evolution of the densities of the atomic conden- 
sate n c and the molecular fraction rib for different evolu- 
tion times of the magnetic field pulse sequence in Fig. Q 
(Evolve = 10, 11, 12, 13 /xs). The initially pure homogeneous 
85 Rb condensate has a density of no = 3.9 x 10 12 cm~ 3 . 

of n^/no = 6%. The small oscillations in n c and n\, 
during the evolution period in Fig. |3| indicate the small 
(n [a(-Bevoive)] 3 = 8 x 10 -5 ) but finite overlap between 
the atomic condensate and the molecular fraction. At 
the lowest value of the field pulses in Fig. ^ one obtains 
no[a(B m i n )] 3 — 0.27 and the observable that determines 
the number of bound pairs in a weakly interacting gas 
even yields rib > fk>/2 (not shown explicitly in Fig. |3J). 
This clearly indicates the significant overlap between the 
molecular wave function of a pair of atoms with its sur- 
rounding gas atoms. The gas is then described by a 
strongly correlated non-stationary many-body state. 

The entire dynamics of the different components of the 
gas can be understood in an intuitive way when the lin- 
ear ramps of the two field pulses in Fig. ^ are idealized 
as sudden switches of the magnetic field strength. The 
first field pulse then shifts the virtually uncorrelated ini- 
tial condensate into a many-body state with the crucial 
binary correlations. The overlap of this many-body state 
with the multi-channel bound state of the pairs of atoms 
at the field strength i3 evo i ve determines the molecular 
fraction in the evolution period. The overlap with the 



corresponding excited binary scattering states provides a 
first burst fraction 0, of correlated pairs of compar- 
atively hot atoms. In the evolution period the atomic 
condensate and the molecular component are virtually 
orthogonal and evolve coherently. At the end of the evo- 
lution period the difference in phase between the ampli- 
tude of the atomic condensate and the molecular fraction 
is Aip = Et,(B cvo i vc )t evo i ve / H. The second field pulse gets 
both components to overlap and, thereby, probes their 
phase difference. The second burst fraction as well as 
the final atomic condensate and molecular component 
thus exhibit an interference depending on Aip. As the 
field pulses are chosen as mirror images in Fig.^ at con- 
structive interference, nb(*an) can be expected to result 
in about twice the molecular density in the evolution pe- 
riod. 

We have shown in this letter that highly excited di- 
atomic bound states produced in recent experiments 
are characterized by a large spatial extent that by far ex- 
ceeds the size of all known ground state molecules [l7|. 
The corresponding 85 Rb2 wave functions are strongly 
dominated by their {—2,-2} open channel component. 
We have analyzed the non-adiabatic association mecha- 
nism of Ref. on the basis of a microscopic quantum 
mechanical many-body description of the gas. Our pre- 
dicted molecular fraction of 6% in the evolution period 
of the pulse sequence Q exceeds the results of Ref. 
by more than an order of magnitude. This provides a 
significantly better perspective for proposed experimen- 
tal schemes to detect the molecules. Corrections imposed 
by the atom trap 7] do not affect the orders of magnitude 
reported here. Their systematic study in connection with 
possible improvements of molecular production schemes 
will be subject to future work. 
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a Royal Society Wolfson Merit Award holder. 



[1] R.H. Wynar et al, Science 287, 1016 (2000). 

[2] J. Stenger et al, Phys. Rev. Lett. 82, 2422 (1999). 

[3] S.L. Cornish et aZ.,Phys. Rev. Lett. 85, 1795 (2000). 

[4] E.A. Donley et al, Nature (London) 417, 529 (2002). 

[5] S.J.J.M.F. Kokkelmans and M.J. Holland, 
Phys. Rev. Lett. 89, 180401 (2002). 

[6] M. Mackie, K.-A. Suominen, and J. Javanainen, 
Phys. Rev. Lett. 89, 180403 (2002). 

[7] T. Kohler, T. Gasenzer, and K. Burnett, Phys. Rev. A 
67, 013601 (2003). 

[8] References 0, separate out physical observables as- 
sociated with diatomic bound states in terms of a phe- 
nomenological quantum field that is included in the 
Hamiltonian. The molecular fraction reported in [f| is 
less than 0.5% throughout the evolution period of the 
pulse sequence and rises steeply at the final ramp. The 



5 



results in this letter will exceed these predictions by more 
than an order of magnitude. 
[9] E. Braaten, H.-W. Hammer, and M. Kusunoki, e-print 
arXive cond-mat/0301489 
[10] Private communication from Eleanor Hodby and Neil 
Claussen. 

[11] T. Kohler and K. Burnett, Phys. Rev. A 65, 033601 
(2002). 

[12] This is the long range asymptotic form of the free Green's 
function evaluated at the near resonant binding energy 
E h = -h 2 /ma 2 (see, e.g., 0). 



[13] J.L. Roberts et al, Phys. Rev. A 64, 024702 (2001). 
[14] F.H. Mies, E. Tiesinga, and P.S. Julienne, Phys. Rev. A 

61, 022721 (2000). 
[15] J. Fricke, Ann. Phys. (N.Y.) 252, 479 (1996). 
[16] For a general description of the detection of composite 

particles see, e.g., J. D. Dollard, J. Math. Phys. 14, 708 

(1973). 

[17] The bond length of the largest known diatomic ground 
state molecule 4 He2 is of the order of 100 aeohr [see 
R.E. Grisenti et al, Phys. Rev. Lett. 85, 2284 (2000)]. 



